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We report on observations of the Large Magellanic Cloud with the Fermi Gamma-Ray Space Telescope. The 
LMC is clearly detected with the Large Area Telescope (LAT) and for the first time the emission is spatially 
well resolved in gamma-rays. Our observations reveal the massive star forming region 30 Doradus as a bright 
source of gamma-ray emission in the LMC. The observations furthermore show that the gamma-ray emission 
correlates little with the gas density of the LMC. Implications of this finding will be discussed. 



1. Introduction 

Since the early days of high-energy gamma-ray as- 
tronomy it has been clear that the gamma-ray flux 
received at Earth is dominated by emission from 
the Galactic disk [1]. This emission is believed to 
arise from cosmic-ray interactions with the interstel- 
lar medium, which at gamma-ray energies ^100 MeV 
are dominated by the decay of tt^ produced in colli- 
sions between cosmic-ray nuclei and the interstellar 
medium |2|]. Further contributions are from cosmic- 
ray electrons undergoing inverse Compton scattering 
off interstellar soft photons and Bremsstrahlung losses 
within the interstellar medium. Gamma-ray observa- 
tions thus have the potential to map cosmic-ray accel- 
eration sites in our Galaxy which may ultimately help 
to identify the sources of cosmic-ray acceleration. 

Nearby galaxies have the advantage of being viewed 
from outside and so line of sight confusion, which com- 
plicates studies of emission from the Galactic disk, is 
diminished. This advantage is however somewhat off- 
set by the limitations by the angular resolution and 
sensitivity of the instrument. The Large Magellanic 
Cloud (LMC) is thus an excellent target for studying 
the link between cosmic-ray acceleration and gamma- 
ray emission since the galaxy is nearby {D ^ 50 kpc) 
d, m, has a large angular extent of ~ 8°, and is seen 
at a small inclination angle of i ^ 20° — 35° [H, @] 
that avoids source confusion. In addition, the LMC is 
relatively active, housing many supernova remnants, 
bubbles and superbubbles, and massive star forming 
regions that are all potential sites of cosmic-ray accel- 
eration 0,11,11. 

The EGRET telescope onboard the Compton 
Gamma-Ray Observatory (CGRO, 1991-2000) was 
the first to detect the LMC p^]. Due to EGRET's lim- 
ited angular resolution and limited sensitivity, details 
of the spatial structure of the gamma-ray emission 
could not be resolved, yet the observations showed 
some evidence for the spatial distribution being con- 
sistent with the morphology of radio emission. The 
Large Area Telescope (LAT) aboard Fermi is provid- 
ing now for the first time the capabilities to go well 
beyond the study of the integrated gamma-ray flux 



from the LMC PJj,ll2[. In this contribution we present 
our first in-depth analysis of the LMC galaxy based 
on 11 months of continuous sky survey observations 
performed with Fermi/LAT. We put a particular em- 
phasis on the determination of the spatial distribution 
of the gamma-ray emission, that, as we will show, re- 
veals the distribution of cosmic rays in the galaxy. A 
more detailed discussion of the observations and their 
implications is given in [13,]. 



2. Observations 
2.1. Data preparation 

The characteristics and performance of the LAT 
aboard Fermi are described in detail by fi^. The 
data used in this work amount to 274.3 days of con- 
tinuous sky survey observations over the period Au- 
gust 8th 2008 - July 9th 2009 during which a total 
exposure of ~ 2.5 x 10^^ cm^ s (at 1 GeV) is ob- 
tained for the LMC. Events satisfying the standard 
low-background event selection ( "Diffuse" events) [l5[ 
and coming from zenith angles < 105° (to greatly re- 
duce the contribution by Earth albedo gamma rays) 
are used. To further reduce the effect of Earth albedo 
backgrounds, the time intervals when the Earth was 
appreciably within the field of view (specifically, when 
the center of the field of view was more than 47° from 
the zenith) are excluded from this analysis. Further- 
more, time intervals when the spacecraft was within 
the South Atlantic Anomaly are also excluded. We 
further restrict the analysis to photon energies above 
200 MeV; below this energy the effective area in the 
"Diffuse class" is relatively small and strongly depen- 
dent on energy. All analysis is performed using the 
LAT Science Tools package, which is available from 
the Fermi Science Support Center, using P6_V3 post- 
launch instrument response functions (IRFs). These 
take into account pile-up and accidental coincidence 
effects in the detector subsystems that are not consid- 
ered in the definition of the pre-launch IRFs. 

At the Galactic latitude of the LMC {b ^ -33°), 
the gamma-ray background is a combination of ex- 
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Figure 1: Background subtracted counts map of the LMC region. The counts map has been adaptively smoothed [T^ 
with a signal-to-noise ratio of 5 in order to reveal significant structures at all possible scales while supressing the noise 



that arises from photon counting statistics. The white line shows the N(Hl) 
to indicate the extent of the gaseous disk. 



10^^ H cm ^ iso column density contour 



tragalactic and Galactic diffuse emissions, some resid- 
ual instrumental background and a number of point 
sources that primarily are associated to blazars. We 
model the background using components for the dif- 
fuse Galactic and the extragalactic and residual in- 
strumental backgrounds and 6 point sources that all 
are associated with known blazars. The Galactic com- 
ponent is based on the LAT standard diffuse back- 
ground model gll_iein_v02 for which we keep the 
overall normalization as a free parameter. The ex- 
tragalactic and residual instrumental backgrounds are 
combined into a single component which has been 
taken as being isotropic. The spectrum of this compo- 
nent is determined by fitting an isotropic component 
together with a model of the Galactic diffuse emis- 
sion and point sources to the data. Also here we leave 
the overall normalization of the component as a free 
parameter. The 6 background blazars are modelled 



as point sources with power-law spectral shapes. The 
flux and spectral power-law index of each source are 
left as free parameters of our background model and 
their values were determined from likelihood analysis. 

2.2. Spatial distribution 

2.2.1. Counts map 

To investigate the spatial distribution of gamma- 
ray emission toward the LMC we show in Fig. [T] an 
adaptively smoothed background subtracted counts 
map of the LAT data. The map clearly shows ex- 
tended emission that is spatially confined to within 
the LMC boundaries which we trace by the iso col- 
umn density contour A^h = 10^^ H cm~^ of neutral 
hydrogen in the LMC [1^]. The total number of ex- 
cess 200 MeV - 20 GeV photons above the background 
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Figure 2: Zoom into a 3° x 3° large region of the 
background subtracted counts map around the central 
star cluster R 136 in 30 Dor. Stars show the locations of 
R 136 and of the pulsars PSR J0540-6919 and 
PSR J0537-6910. The circles show the 95% containment 
radius of sources SI and S3 of the point source model, 
and of source G2 of the 2D Gaussian shaped model (see 
text). During month 4 of our dataset, a flaring point 
source occured near 30 Dor, and we indicate the 95% 
containment radius of this source; the diamond shows the 
location of the possible counterpart RX J0546.8— 6851 of 
the flaring source. 

in the LMC area amounts to ~ 1550 counts whereas 
the background in the same area amounts to ~ 2440 
counts. With these statistics, the extended gamma- 
ray emission from the LMC can be resolved into 
several components. The brightest emission feature 
is located near (aj2ooo, ^J200o) ~ (05^0"^, -69°150, 
which is close to the massive star-forming region 
30 Doradus (30 Dor) that houses the two Crab- 
like pulsars PSR J0540-6919 and PSR J0537-6910 
[13, 13 • Excess gamma-ray emission is also seen to- 
ward the north and the west of 30 Dor. These bright 
regions are embedded into a more extended and dif- 
fuse glow that covers an area of approximately 5° x 5°. 
Figure [2] present a zoom of Fig. [1] of the 30 Dor region 
over which we overlay potential sources of gamma-ray 
emission. 

2.2.2. Model fitting 

As next step, we assess the spatial distribution of 
the LMC emission using (1) simple parametrized geo- 
metrical models of the gamma-ray intensity distribu- 
tion, and (2) spatial templates that trace the interstel- 



lar matter distribution in the LMC. We assume power- 
law spectral distributions for all models and keep the 
total flux and power law index as free parameters. 
We adjust the spatial and spectral parameters of the 
models using a binned maximum likelihood analysis 
with spatial bins of 0.1° x 0.1° and 60 logarithmically 
spaced energy bins covering the energy range 200 MeV 
- 20 GeV. We quantify the goodness-of-fit using the 
so-called Test Statistic (TS) which is defined as twice 
the difference between the log-likelihood £i that is ob- 
tained by fitting the model on top of the background 
model to the data, and the log-likelihood Co that is 
obtained by fitting the background model only, i.e. 

TS = 2(/:i-/:o). 

First, we examine if the gamma-ray emission from 
the LMC can be explained with a combination of in- 
dividual point sources. For this purpose we add suc- 
cessive point sources to our model and optimize their 
locations, fluxes and spectral indices by maximizing 
the likelihood of the model. We stop this procedure 
once the TS improvement after adding a further point 
source drops below 25. This happens after we added 
5 point source to our model, resulting in TS = 1089.3. 
The left panel of Fig.[3]shows the corresponding model 
counts map. 

Second, instead of using point sources we repeat 
the procedure with 2D Gaussian shaped intensity pro- 
files to build a geometrical model that is more appro- 
priate for extended and diffuse emission structures. 
We again stop the successive addition of 2D Gaussian 
shaped sources once the TS improvement after adding 
a further source drops below 25. This occurs after two 
2D Gaussian shaped sources have been added to the 
model, resulting in TS = 1122.6. We show the cor- 
responding model counts map in the right panel of 
Fig. [3l Obviously, the 2D Gaussian model provides 
a larger TS than the point source model, suggesting 
that it better fits the data despite the smaller num- 
ber of free parameter (10 for the 2D Gaussian model 
compared to 20 for the point source model). It is thus 
more likely that the LMC emission is indeed diffuse 
in nature, or alternatively, composed of a large num- 
ber of unresolved and faint sources that can not be 
detected individually by Fermi/ L AT. 

Third, we fit the LAT data to gas maps of neu- 
tral atomic hydrogen (Hi) pj], molecular hydrogen 
(H2) 0, and ionized hydrogen (Hii) 0. Fitting 
the Hi and H2 maps result in TS values of 771.8 and 
824.3, respectively, that are considerably worse than 
those obtained for the geometrical models. Appar- 
ently, the H I and H2 maps provide rather poor fits to 
the data, indicating that the distribution of gamma 
rays does not follow the distribution of neutral hy- 
drogen in the LMC. This can already be seen from 
the left and mid panels of Fig. [4] which show the gas 
maps after convolution with the LAT point spread 
function. While the gamma-ray emissivity is high- 
est in 30 Dor and the northern part of the galaxy. 
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Figure 3: Best fitting point source model (left panel) and 2D Gaussian shape model (right panel) that has been 
convolved with the LAT point spread function. The point source model consists of 5 point sources (S1-S5) while the 2D 
Gaussian shaped model consists of a broad component covering a large fraction of the LMC disk (Gl) and a narrow 
component near 30 Dor (G2). 




Figure 4: Gas tracer maps convolved with the LAT point spread function and scaled using maximum likelihood model 
fitting to the data. Panels show Hi (left panel), H2 (mid panel), and Hll (right panel). 



the gas maps show a bright ridge that runs over 
~ 3° along aj2ooo ^ 05^40"^ which coincides with the 
most prominent region of ^^CO emission tracing gi- 
ant molecular clouds in the LMC [21]. Roughly 20% 
of the total gas mass in the LMC is confined into this 
ridge [22], yet comparison with Fig. [T] shows that most 
of the ridge is not luminuous in high-energy gamma 
rays. 



The H II map, on the other hand, gives TS = 1110.1, 
which is very close to the TS values of the geometrical 
models. The H 11 map thus provides the best fit among 
all of the gas maps to the LAT data, and the right 
panel of Fig. [4] indeed shows that the model is very 
similar to the 2D Gaussian shaped model and also fol- 
lows closely the observed distribution of gamma rays 
(Fig.HD. 



2.3. Emissivity spectrum 

To determine the spectrum of the gamma-ray emis- 
sion from the LMC independently from any assump- 
tion on the spectral shape, we fit our data in 6 log- 
arithmically spaced energy bins covering the energy 
range 200 MeV - 20 GeV. We obtain the total spec- 
trum of the LMC by fitting the Hii template to the 
data. We also obtain separate spectra for the LMC 
disk (Gl) and for 30 Dor (G2) by fitting the 2D Gaus- 
sian shaped model to the data. 

To determine the integrated gamma-ray flux we fit 
exponentially cut off power law spectral models of the 
form N{E) = k (E/Eo)-^ exp{-E/Ec) to the data. 
We make these fits by means of a binned maximum 
likelihood analysis over the energy range 200 MeV - 
20 GeV. This results in an extrapolated > 100 MeV 
photons fiux of (2.6 ±0.2) x 10"'^ ph cm ^ s ^ for the 
Hii template which corresponds to an energy flux of 
(1.6 ± 0.1) X 10~^^ erg cm~^ s~^ (systematic uncer- 
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Figure 5: Differential average gamma-ray emissivity for 
the total LMC emission (black dots) and for the LMC 
disk only (red dots) compared to that of the Galactic 
local interstellar medium determined by [23] (grey data 
points). The solid line shows the predicted gamma-ray 
emissivity computed in the framework of a one-zone 
model for the LMC disk O. The other lines show the 
contributions of 7r°-decay (long dashed), Bremsstrahlung 
(short dashed), and inverse Compton emission (dotted). 



tainties in these estimates amount to less than 16%). 

Dividing our spectra by the spatially integrated hy- 
drogen column density J = (3.6±1.2) x 10^^ H- 
atom cm~^ sr of the LMC provides us with the differ- 
ential gamma-ray emissivity per hydrogen atom. We 
compare the resulting emissivity for the total LMC 
emission and the LMC disk only in Fig. [5] to the one 
obtained for the Galactic local interstellar medium by 
(23| . Apparently, the local Galactic emissivity is be- 
tween ~ 2 to ^ 4 times larger than the average emis- 
sivity of the LMC. 

We compare the differential gamma-ray emissivi- 
ties to a one-zone model of cosmic-ray interactions 
with the interstellar medium that takes into ac- 
count TT^ decay following proton-proton interactions, 
Bremsstrahlung from cosmic-ray electrons and inverse 
Compton scattering of cosmic-ray electrons on LMC 
optical and infrared photons and cosmic microwave 
background photons [l^]. Fitting this spectral model 
to our data using a binned maximum likelihood anal- 
ysis gives an average cosmic-ray enhancement fac- 
tor of Tc = 0.31 ± 0.01 for the entire LMC, and of 
Tc = 0.21 ± 0.01 for the LMC disk only. System- 
atic errors due to uncertainties in the effective area 
of the instrument amount to ±0.02. An additional 
systematic error of —23% to +42% comes from the 
uncertainty in the total gas mass of the LMC, which 
largely dominates the statistical and systematic mea- 
surements errors. 



3. Cosmic-ray density distribution 

To reveal the sites of cosmic-ray acceleration in the 
LMC we map the cosmic-ray density variations in 
the galaxy by computing the gamma-ray emissivity 
as function of position. We do this by dividing our 
background subtracted counts map by the A/'(H) map 
after convolution of the latter with the LAT instru- 
mental response function. We normalise A^(H) to a 
total LMC hydrogen mass of 7.2 x 10^ Mq that takes 
into account the possible presence of dark gas that 
is not seen in radio surveys of Hi [g^. We adap- 
tively smooth [14] the counts maps and used the re- 
sulting smoothing kernel distribution to smooth also 
the convolved A^(H) map before the division to reveal 
significant structures at all possible scales, while su- 
pressing the noise that arises from the limited photon 
counting statistics. The resulting emissivity maps are 
shown in Fig. [6l We superimpose on the images the 
interstellar gas distribution, as traced by A^(H), con- 
volved with the LAT instrumental response function, 
and also show the locations of potential particle ac- 
celeration sites, such as pulsars, supernova remnants, 
Wolf-Rayet stars and supergiant shells. 

Figure [6] reveals that the cosmic-ray density varies 
considerably over the disk of the LMC. The gamma- 
ray emissivity is highest in 30 Dor and the northern 
part of the galaxy, while the southern part and in par- 
ticular the dense ridge of gas south of 30 Dor seems 
basically devoid of cosmic rays. These large variations 
confirm our earlier findings that the gamma-ray emis- 
sion correlates little with the gas density in the LMC. 
Figure [6] suggests further that the cosmic-ray den- 
sity correlates with massive star forming tracers, and 
in particular Wolf-Rayet stars and supergiant shells. 
This finding is corroborated by the good fit of the H ii 
gas map, which is probably the most direct tracer of 
massive star forming regions within a galaxy. 

Thus, the gamma-ray emissivity maps of the LMC 
support the idea that cosmic rays are accelerated in 
massive star forming regions as a result of the large 
amounts of kinetic energy that are input by the stel- 
lar winds and supernova explosions of massive stars 
into the interstellar medium. Our data reveal a rela- 
tively tight confinement of the gamma-ray emission to 
star forming regions, which suggests a relatively short 
diffusion length for GeV protons. 



4. Conclusions 

Observations of the LMC by Fermi/LAT have for 
the first time provided a detailed map of high-energy 
gamma-ray emission from that galaxy. Our analysis 
reveals the massive star forming region 30 Doradus as 
bright source of gamma-ray emission in the LMC in 
addition to fainter emission regions found in the north- 
ern part of the galaxy. The gamma-ray emission from 
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Figure 6: Integrated > 100 MeV emissivity maps of the LMC in units of 10 ph s ^ sr ^ H-atom ^. An adaptive 
smoothing with a signal-to-noise ratio of 5 has been apphed to reduce statistical fluctuations. The white line shows the 
N(Hl) = 10^^ H cm ^ iso column density contour to indicate the extent of the gaseous disk. Symbols indicate the 
locations of pulsars (pluses) from ATNF catalogue version 1.36 [24], supernova remnants (diamonds) from Rosa 
Williams web page http://www.astro.illinois.edu/projects/atlas/index.html, Wolf-Rayet stars (stars) from the fourth 
catalogue of and supergiant shells (circles) from [2^. 



the LMC shows very little correlation with gas density. 
A much better correlation is seen between gamma-ray 
emission and massive star forming regions, as traced 
by the ionizing gas, Wolf-Rayet stars and supergiant 
shells, and we take this as evidence for cosmic-ray ac- 
celeration in these regions. This correlation supports 
the idea that cosmic rays are accelerated in massive 
star forming regions as a result of the large amounts 
of kinetic energy that are input by the stellar winds 
and supernova explosions of massive stars into the in- 
terstellar medium. 

Continuing observations of the LMC with 
Fermi/LAT in the upcoming years will provide 
the photon statistics to learn more about the origin 
of the gamma-ray emission from that galaxy. Better 
statistics will help in identifying more individual 



emission components and may help to separate true 
point sources from the more diffuse emission that we 
expect from cosmic-ray interactions. 
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